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EXECUTIVE  SUMMARY 


PROBLEM 

The  success  of  accurately  tracking  an  aerospace  object  is  strongly  dependent  on  how 
well  the  position,  vetocity,  and  acceleration  of  the  object  is  known.  Understanding 
what  contributes  to  errors  in  these  quantities,  and  how,  is  important  in  predicting 
the  performance  of  a  tracking  system,  or  as  in  the  case  of  a  Single  Integrated  Air 
Picture  (S!AP),  a  distributed  tracking  system,  A  previous  JSSEO  Technical  Report  [1] 
investigated  position  estimate  errors  as  a  function  of  data  registration  and  dynamic  errors. 
This  Technical  Report  is  a  complement  to  that  earlier  analysis  and  outlines  a  similar 
methodology  for  understanding  and  investigating  track  velocity  estimate  errors.  The 
study  of  velocity  estimate  errors  is  relevant  not  only  to  the  general  problem  of  aerospace 
object  tracking,  but  also  to  the  analysis  of  specific  operational  threads.  For  example,  in  an 
Integrated  Fire  Control  (IFC)  engagement,  the  interceptor  will  use  positional  and  velocity 
information  about  the  target.  Therefore,  knowledge  of  the  errors  in  both  the  position  and 
velocity  estimates  is  important  to  predicting  system  performance. 

OBJECTIVES 

This  Technical  Report  presents  an  analytical  framework  for  quantifying  track  velocity 
errors  as  a  function  of  data  registration  and  dynamic  errors,  The  methodology  is 
general  and  can  be  applied  to  fand-based,  sea-based,  and  airborne  sensor  systems. 
The  technique  provides  a  straightforward  means  of  identifying  the  major  elements  that 
contribute  to  error  in  track  velocity  estimates  and  in  quantitying  their  contribution  to  track 
velocity  error. 

APPROACH 

The  basis  of  the  analysis  is  to  develop  the  expression  that  relates  the  velocity  of  the 
aerospace  object  expressed  in  a  sensor  frame  to  the  velocity  of  the  object  relative 
to  the  Earth's  center  expressed  in  the  WGS  84  coordinate  frame,  and  to  consider 
how  measurement,  navigation,  alignment,  tracker,  and  liming  {e.g..  time  stamp,  time 
latency,  time  synchronization,  etc.)  errors  affect  the  estimate  of  track  velocity.  A  linear 
error  analysis  is  used  to  characterize  and  understand  the  track  velocity  estimation 
performance  of  an  arbitrary  sensor  system.  The  process  involves  identifying  all  elements 
that  introduce  error,  tracing  the  propagation  of  these  errors  throughout  the  system,  and 
quantifying  their  effect  on  system  performance.  To  make  the  analysis  more  tractable, 
certain  reasonable  assumptions  are  made.  Finally,  a  nominal  tracking  scenario  is  used 
to  illustrate  the  magnitude  of  the  track  velocity  estimate  error  and  identify  which  system 
errors  are  the  major  contributors  to  the  total  track  velocity  estimate  error, 

FINDINGS 

The  methodology  previously  used  to  quantify  track  position  estimate  error  can  be 
straightfonwardly  extended  to  the  investigation  of  track  velocity  estimate  error.  The 
resulting  track  velocity  error  sensitivities  are  very  similar  in  form  to  their  track  position 
error  counterparts  but  contain  some  additional  complexity.  Comparing  the  results  for 
velocity  estimate  error  to  those  for  position  estimate  error,  we  identify  new  contributors 
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for  (rack  velocity  error:  vehicle  velocity  errors  and  sensor  rate  measurement  errors, 
l0.,  range  rate,  bearing  rate,  and  elevation  rate  (if  they  are  available).  Relative  to  track 
position  error,  track  velocity  error  is  much  less  sensitive  to  the  individual  system  errors. 

CONCLUSIONS 

This  Technical  Report  provides  the  details  of  an  analytical  technique  to  understand  and 
quantify  how  errors  in  t^avigation,  sensor,  timing,  and  tracker  affect  distributed  system 
performance.  A  key  result  of  the  analysis  is  the  determination  of  the  sensitivities  of 
track  velocity  accuracy  to  navigation,  sensor,  tracker,  and  timing  errors.  A  nominal 
tracking  scenario  is  considered  and  both  the  velocity  error  sensitivities  and  the  position 
error  sensitivities  are  computed  and  examined.  In  general,  track  position  error  is  much 
more  sensitive  to  the  various  system  error  terms.  For  the  scenario  considered,  the 
track  velocity  error  sensitivities  to  tracker  acceleration  errors  and  timing  errors  are  the 
greatest.  For  position  error,  the  sensitivities  to  angular-type  errors  and  timing  errors  are 
dominant.  Some  similarities  exist  between  the  velocity  error  sensitivities  and  the  position 
error  sensitivities.  Specifically,  the  position  error  sensitivities  to  navigation  position  errors 
and  sensor-measured  range  errors  are  equivalent  to  the  velocity  error  sensitivities  to 
navigation  velocity  errors  and  sensor-measured  range  rate  errors,  which  are  just  the  time 
rates  of  change  of  navigation  position  and  range. 

Using  the  velocity  error  sensitivities  with  a  nominal  system  error  budget  for  the  chosen 
scenario  yields  a  sample  track  velocity  accuracy  error  budget.  For  this  particular 
scenario,  dynamic  errors  dominate  the  total  track  velocity  error.  These  dynamic  errors 
are  composed  of  tracker  and  timing  error  contributions.  Specifically,  the  tracker  velocity 
and  acceleration  error  contributions  are  the  primary  cause  for  the  large  tracker  error 
contribution. 

We  will  to  incorporate  the  results  of  this  analysis  into  the  JSSEO  Parametric  Track 
Accuracy  Model.  This  will  increase  the  scope  of  the  model  and  allow  for  additional 
studies  and  analysis  to  be  conducted.  For  example,  the  results  of  the  parametric 
model  and  the  system  specific  error  parameters  can  be  applied  to  using  and  testing  the 
Integrated  Architecture  Behavior  Model  (lABM)  within  the  Joint  Distributed  Engineering 
Plant  (JDEP)  Technical  Framework.  In  particular,  the  system  specific  error  parameters 
define  the  statistical  error  distribution  for  each  error  term.  The  parameter  values  coupled 
with  the  corresponding  distribution  define  the  level  of  bias  that  should  be  introduced  into 
every  sensor  that  is  associated  with  an  instantiation  of  the  lABM.  The  magnitude  of  the 
contribution  of  each  error  to  track  velocity  error  provides  a  hierarchy  for  the  sequence 
of  adding  errors  and  systems  into  the  lABM  for  testing.  Since  larger  error  contributors 
dominate  smaller  contributors,  the  errors  with  smaller  contributions  should  be  introduced 
into  the  lABM  first  when  only  one  system  is  being  considered.  A  similar  hierarchy  should 
also  be  applied  to  testing  the  iABM  within  a  distributed  system.  The  more  accurate 
systems  should  be  introduced  into  the  distributed  system  before  the  less  accurate 
systems.  The  results  of  the  error  model  analysis  also  provide  examples  of  track  velocity 
error  that  the  IABM  can  be  checked  against. 
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RECOMMENDATIONS 

The  results  indicate  that  track  velocity  error  contributors  are  less  important  than  track 
position  error  contributors  in  meeting  distributed  system  capability  requirements. 
Therefore,  the  lABM  development  activities  should  focus  on  capturing  the  position  error 
contributors  for  the  Configuration  05  build. 
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1  INTRODUCTION 

1.1  Problem  Statement 

[n  a  previous  report  [1],  a  methodology  was  developed  to  investigate  and  identify  the 
major  data  registration  and  timing  elements  that  contribute  to  error  in  track  positional 
estimates.  The  analysis  contained  in  that  report  provides  the  basis  for  the  ongoing 
JSSEO  Parametric  Track  Accuracy  Model  development.  This  Technical  Report  extends 
that  analysis  to  include  the  impact  of  system  data  registration,  timing,  and  tracker  errors 
on  track  veiocity  estimates.  Analogous  to  the  investigation  of  track  position  error,  the 
methodology  provides  useful  Insight  into  how  the  various  elements  contribute  to  track 
velocity  estimates.  In  addition,  the  resuits  from  this  analysis  of  track  velocity  error  can  be 
incorporated  into  the  Parametric  Track  Accuracy  Model  thereby  enhancing  its  scope. 

1 .2  Motivation 

The  exchange  of  track  positional  and  velocity  information  is  critical  to  successfully 
integrating  air  pictures  from  various  tracking  systems.  Of  equal  importance  is  the 
exchange  of  the  uncertainties  in  those  track  position  and  velocity  estimates.  For  a 
peer-to-peer  environment,  the  exchange  of  track  data  and  its  associated  uncertainty 
between  the  participants  will  be  critical  to  the  distributed  system  performance.  It  is  also 
important  to  understand  the  origin  and  behavior  of  the  uncertainty  such  that  system 
performance  can  be  analyzed  and  improvements  identified. 

In  coordinating  the  participating  systems  in  a  network-centric  environment  such  as 
Unk-16,  a  Track  Quality  (TQ)  Number  is  used  to  quantify  track  position  uncertainty 
by  reducing  a  three-dimensional  position  error  covariance  to  a  single  number  that  is 
transmitted  along  with  the  track  data  [2],  Currently,  only  horizontal  position  errors 
enter  into  the  calculation  of  TQ,  however,  future  improvements  may  lead  to  additional 
information  contributing  to  the  TQ  Number.  For  instance,  the  tracking  algorithms  used 
by  the  participating  systems  will  often  include  the  three  components  of  velocity  in  their 
state  vectors.  Covariance  information  for  the  velocity  error  could  therefore  be  very  easily 
included  in  the  calculation  of  TQ.  An  understanding  of  what  contributes  to  track  velocity 
error  is  necessary  to  determine  how  Track  Quality  would  be  affected  by  the  addition  of 
velocity  error  information.  Such  knowledge  also  enables  investigations  into  other  relevant 
areas,  e.g,,  tracking  algorithms. 

The  study  of  velocity  estimate  errors  is  relevarrt  not  only  to  the  general  problem  of 
aerospace  object  tracking,  but  also  to  the  analysis  of  specific  operational  threads. 

For  example,  in  an  Integrated  Fire  Control  (IFC)  engagement,  the  interceptor  will  use 
positional  and  velocity  information  about  the  target.  Therefore,  knowledge  of  the  errors  in 
both  the  position  and  velocity  estimates  is  important  to  predicting  system  performance. 
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1.3  Technical  Approach 


This  report  largely  follows  the  same  structure  as  the  earlier  report  on  track  position  error 
analysis  [1],  Therefore,  some  famiiiarity  with  the  content  of  that  report  is  assumed  to 
reduce  redundancy  here.  However,  Appendix  A  contains  the  position  error  analysis 
results  necessary  for  the  present  treatment  of  velocity  error.  The  basis  of  the  analysis  is 
to  develop  the  expression  that  relates  the  velocity  of  the  aerospace  object  expressed  in  a 
sensor  frame  to  the  velocity  of  the  object  relative  to  the  Earth’s  center  expressed  in  the 
WGS  84  coordinate  frame,  and  to  consider  how  measurement,  navigation,  alignment, 
timing  (e.g.,  time  stamp,  time  latency,  time  synchronization,  etc.),  and  tracker  errors  affect 
the  estimate  of  track  velocity. 

The  analysis  is  general  and  can  be  applied  to  land-based,  sea-based,  and  airborne 
sensor  systems.  As  with  studying  the  position  error,  a  linear  error  analysis  is  used  to 
characterize  and  understand  the  track  velocity  estimation  performance  of  an  arbitrary 
sensor  system.  The  process  involves  idenli^ing  ait  elements  that  introduce  error, 
tracing  the  propagation  of  these  errors  throughout  the  system,  and  quantifying  their 
effect  on  system  performance.  In  the  analysis,  the  system  errors  are  assumed  to  be 
normal  random  variables  represented  by  a  deterministic  mean  and  a  variance.  This 
type  of  analysis  has  many  benefits,  some  of  which  include:  quantifying  how  individual 
subsystem  (e.g.,  navigation)  as  well  as  individual  elements  (e.g.,  latitude  error)  contribute 
to  overall  system  performance,  identifying  potential  improvement  areas,  predicting 
system  performance,  and  enabling  trade-off  studies. 
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2  VELOCITY  ERROR  ANALYSIS  FOR  AEROSPACE  OBJECTS 


2.1  The  Aerospace  Object  Velocity  Vector  Processing  Chain 

In  the  following  series  of  steps,  the  nominal  (error-free)  aerospace  object  velocity  vector 
processing  chain  Is  presented.  The  starting  point  of  the  processing  chain  is  the  sensor 
measurement  of  the  object:  range  p,  bearing  ii>«,  and  elevation  6^.  The  ending  point  will 
be  the  object  velocity  vector  expressed  in  Earth  coordinates.  The  formulation  follows 
directly  from  the  position  analysis  (see  Appendix  A). 


Step  1 :  The  sensor  system  measures  the  sensor-object  displacement  vector  in  sensor 
coordinates. 

■  ■ 


=  P 


(2.1) 


where  the  notation  {  la  or  [  ],  represents  the  vector  expressed  in  the  .i'-coordlnate  frame 
and  C„  and  S„  are  shorthand  for  coso  and  aiiuj.  Differentiating  (Ar^t)*  with  respect  to 
time  yields  the  time  rate  of  change  of  the  sensor-object  displacement  vector  (in  the 
sensor  frame), 


pS,„C\  - 

-  f)Se,  -  0X'b, 


(2.2) 


The  term  p  is  the  range  rale,  the  term  is  the  elevation  time  rate  of  change,  and  the 
term  is  the  bearing  time  rate  of  change.  lA/hile  a  sensor  may  provide  a  direct  measure 
of  the  range  rate,  direct  sensor  measurements  of  the  angular  rates  of  change  are  very 
unlikely.  Therefore,  if  we  only  consider  the  role  of  range  rate,  (Afjf).,  reduces  to 


(2,3) 


which  mimics  the  expression  above  for  ( Af,,)^. 


Step  2;  For  position,  the  sensor-object  displacement  vector  Is  transformed  to  body 
coordinates  using  the  sensor-body  rotation  matrix,  Rbh, 


{Ar„)j,  =  /?*.,(  AT,, )„,  (2.4) 

Taking  the  time  derivative  of  this  expression  yields  the  time  rate  of  change  of  the 
sensor-object  vector  (in  the  body  frame): 

{ Ar „ )ft  -  T4s  I A'r,{ ).,  +  ( Ar,: )« •  (2.5) 


The  term  iS  known  from  the  position  analysis  (Equation  A, 5),  and  the  terms  (A'rj,), 
and  tArsi)„  are  given  in  Step  1  (Equations  2.1  and  2,2),  The  new  term  for  the  velocity 
analysis  is  which  is  the  time  rate  of  change  of  the  sensor-body  rotation  matrix.  If  we 
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make  the  simplifying  assumption  that  the  sensor  frame  is  not  rotating  relative  to  the  body 
frame,  then  =  0  and  the  above  expression  reduces  to 

(2*6) 

This  assumption  is  expected  to  be  valid  in  general,  even  for  systems  such  as  rotating 
radars.  The  reason  is  that  even  though  a  radar  may  be  mechanically  rotating  as  part  of 
its  measurement  operation,  Its  measurement  frame  itseif  is  stationary. 

Step  3;  For  position,  the  body-sensor  iever  arm  vector,  is  added  to  yield  the 

body-object  displacement  vector  in  body  coordinates, 

(AfTj)ft  =  (A/Tj);,  +  (Arj()t.  (2,7) 

The  time  derivative  of  this  expression  is  simply, 

{ Afti  )fr  =  ( Affc,  )i  +  ( {2.8) 

If  we  make  the  assumption  that  the  sensor  is  rigidly  fixed  to  the  vehicle  body  and  not 
translating  in  any  way.  (Ar),.,)*  =  0,  and 

{ArT,)t  =  (Af,t}i.  (2.9) 


Step  4:  For  position,  the  body— object  displacement  vector  is  transformed  to  local-level 
coordinates  using  the  local-level-body  rotation  matrix,  ffy, 

{A4);  =  RliAr,,),.  (2.10) 

We  take  the  lime  derivative  of  this  expression  to  obtain  the  time  rate  of  change  of  the 
body-object  vector  (in  the  local-levet  frame): 

t  AT,; );  =  ( A,^  k  +  (2.11) 

The  terms  /fjj,'  and  (AFwJi  are  known  from  the  position  analysis  (Equations  A.9  and  A.7), 
and  the  term  _( Arf,,,)^  is  given  by  Equation  2.8  in  Step  3.  The  new  term  for  the  velocity 
analysis  is  which  is  the  time  rate  of  change  of  the  focal-level-body  rotation  matrix. 

If  we  make  the  simplifying  assumption  that  the  vehicle  body  is  not  pitching,  rolling,  or 
yawing  (or  that  the  rates  of  these  angular  motions  are  sufficiently  small),  then  7?^  -r  u  and 
the  above  expression  simplifies  to 

(AFii),  -  (2.12) 

Step  5:  For  position,  the  body-object  displacement  vector  is  transformed  to  Earth 
coordinates  using  the  Earth-local-levef  rotation  matrix,  7?,^. 

(A4}.= (2.13) 
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We  take  the  time  derivative  of  this  expression  to  obtain  the  time  rate  of  change  of  the 
body-object  vector  (in  the  Earth  frame): 

( h  Ri (i 4 )i .  (2. 1 4} 

The  terms  Rj^,  and  (tirhi);  are  known  from  the  position  analysis  (Equations  A.t  1  and 
A.8),  and  the  term  is  given  by  Equation  2,1 1  in  Step  4,  The  new  term  for  the 

velocity  anaiysis  is  which  is  the  time  rate  of  change  of  the  Earth-focal-level  rotation 
matrix.  Here,  we  will  not  assume  is  zero  since  it  is  a  function  of  vehicle  velocity.  The 
Earth-iocal-levet  rotation  matrix  is 


-S^Cx  -.%Sx 

-Sx  Cx 

0 

-C^Cx  -C\,Rx 

“ 

The  time  derivative  of  Ri^  is 


Rir 


~  CxCj  -S^C,X  - 
-  -S,\ 

C^Sx^  +  CxS^^  -C,,Cx^  +  SxSj^ 


tJ 


(2,15) 


(2.16) 


Entering  into  fl,,  are  the  time  derivatives  of  vehicle  latitude  and  longitude,  0  and  A,  which 
are  refated  to  the  north  and  east  velocities  of  the  vehicle,  If  the  vehicle  is  stationary, 

Ru  =  0. 


Step  6:  For  position,  the  navigation  system  position  vector,  is  added  to  yield  the 
object  position  vector  in  WGS  84  coordinates, 


(n).  =  (rj).  +  (2.17) 

The  time  derivative  of  this  expression  is  simply, 

(r'()e  =  (rt)c  -|-  (Ar'bi)s..  (2.18) 

The  term  (Ar*,,)^  is  given  by  Equation  2.14  in  Step  5,  leaving  (7;,)^,  the  time  rate  of 
change  of  the  navigation  system  position  vector,  as  the  new  term.  From  the  positron 
analysis  (Equation  A.1), 


(/W  +  /0C'pC.\ 
[{1.  -  -I-  /h] 


where 

'Vj,I  = 

The  time  derivative  of  [1%]^  takes  the  form: 


■  -fr.^  +  h}{C^Sx^  +  Ca6>)  +  (4„  +  h.]C^Cx  1 

+  /ijiCoCAA  -  SxSo0) 

■1"  (  "1” 

[(1  -  f"  IvVu  +  !l] 

(i 

(2.19) 

(2.20) 

(2.21) 
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where 


The  equality, 

can  be  used  to  simplify 


(R,}r  =  (f'-M  + 

'  -S^Cx  ' 

-h  {  4*  h 

R)  J 
>■ 

+/p 

1 - 

1 _ 

iJ 

L  '^4  J 

where  the  north-south  radius  of  curvature  is 


(2.22) 

(2.23) 


(2.24) 


r„.,  =  r..{l-^")(l  -  (2-25) 

This  expression  provides  the  vehicle  velocity  vector  expressed  in  Earth  coordinates  as  a 
function  of  the  vehicle  geodetic  position  and  velocity.  Notice  that  the  expression  for  the 
vehicle  velocity  is  of  the  same  form  as  the  expression  for  the  error  in  the  position  of  the 
vehicle  (Equation  A, 33}. 


2.2  The  Aerospace  Object  Error  Chain 


The  same  steps  as  in  the  preceding  section  are  used  to  derive  the  measured  quantities 
from  the  fundamental  error-free  quantities.  Measured  quantities  will  be  denoted  with 
a  ~  and  error  quantities  will  be  denoted  with  a  preceding  <5. 

Step  1:  Starting  again  with  the  position  analysis  (Equation  A, 13},  the  error  in  the 
sensor-object  displacement  vector  in  sensor  coordinates  is. 

(Ar^t),,  =  (Arjt)4  +  (^rsij.-i-  (2,26) 

Differentiating  with  respect  to  time  yields  the  time  rate  of  change  of  the  error  in 

the  sensor-object  displacement  vector  (in  the  sensor  frame), 


(flV’ r.)*  =  (Ah,f).,  -  (irVf).,,  (2.27) 


Following  the  position  error  analysis,  p  =  />  i  =  fi  -+  and  tt'.,  -  t’s  -l-  tii'),;-  This 
yields  in  terms  of  the  sensor  measurement  errors: 


(4r^f),  =  ip 

+  dO^ 

~  P^On  i'r.H,  ~  P  ^  ->  Ti'-'fl  0.,  1.>JI  ^0, 

1 

O 

LA 

1 

—pC'o^  + 

■  +  pfKs,„s,,  ■ 

P>C„X\: 

+  f)p 

(.1 

L  j 

pxW,, 

■S'.,.5V 

+  piVs: 

a/v 

1 

■P 

i _ 

It 
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Note,  the  last  three  terms  in  the  above  expression  mimic  the  result  for  (if,*).,  in  the 
position  analysis  (Equation  A.  14).  If  we  orsly  consider  the  range  rate  p  and  the  range  rate 
error  Sp,  the  above  expression  reduces  to 


4 

'■O 

1 

1 _ 

■  ■ 

(4rj,)^  =  -pdOi 

+  p5i;^ 

-T  Spl 

tt 

■ - [ 

1 

which  similarly  resembles  (Equation  A.  14). 

Step  2:  From  the  position  analysis  (Equation  A. 15),  the  error  in  the  sensor— obfect 
displacement  vector  in  body  coordinates  is 


(2,30) 

Differentiating  (ifst)*,  with  respect  to  lime  yields  the  time  rate  of  change  of  the  error  in 
the  sensor-object  displacement  vector  (in  the  body  frame), 

-T  (A?,*)(,  -  {Af,()t.  (2.31) 

The  full  expression  for  (rtf  *)(,  is 

(-ff« )t  =  ).  -F  m. { Af  f [irT, ),.  (2.32) 

Note,  the  first  two  terms  in  the  above  expression  mimic  the  result  for  (^ir^*)*,  in  the  position 
analysis  (Equation  A.16).  If  we  invoke  the  earlier  assumption  that  Rf^  =  0  and  also  that 
=  0,  then  this  expression  simplifies  to 

=  RiMKi),  +  ^i.,(  Af,,).,.  (2.33) 

The  terms  T4.,  and  are  provided  by  the  position  analysis  (Equations  A.5  and  A.17), 
(4f,,)„  is  given  by  Equation  2.28  in  Step  1,  and  {A7~a()j,  is  provided  by  Equation  2.2  in 
Step  1  of  Section  2.1, 

Step  3:  From  the  position  analysis  (Equation  A,  19),  the  error  in  the  body-object 
displacement  vector  in  body  coordinates  is 

(Afiij)i,  =  (Af,,)*,  -f-  (t5f{,()£,.  (2.34) 

Differentiating  (A7f,,)i  with  respect  to  time  yields  the  time  rate  of  change  of  the  error  in 
the  body-object  displacement  vector  (in  the  body  frame), 

)h  =  { Afi„ -  ( Af'o*  )fr.  (2.35) 

Since  we  are  assuming  that  there  is  no  time  rate  of  change  of  the  body-sensor  lever 
arm,  i.e.,  {At%)(,  =  f),  we  will  also  assume  (i'fiJti  =  d-  9i\^ing 

{af(„)fr-{4f„,),..  (2,36) 
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with  being  provided  by  Equation  2.32  in  Step  2. 

Step  4:  From  the  position  analysis  (Equation  A.21),  the  error  in  the  body-object 
displacement  vector  in  local-level  coordinates  is 

=  (Ar«),  +  (2,37) 

Differentiating  (AFk);  with  respect  to  time  yields  the  time  rate  of  change  of  the  error  in 
the  body-object  displacement  vector  (in  the  local-level  frame), 

(tfFii),  =  (An,),  -  (AF,,,),.  (2.33) 

The  full  expression  for  (<irw.),  takes  the  form 

Ar-i,)i  +  4(;(Aft,  k  +  Rlmik-  (2-39) 

Note,  the  first  two  terms  in  the  above  expression  mimic  the  result  for  (jSfftt),  in  the  position 
analysis  (Equation  A. 22],  If  we  invoke  the  earlier  assumption  that  =  0  and  also  that 
=  0,  then  this  expression  simplifies  to 

),  =  )i  4  ( A Fw )t.  (2.40) 

The  terms  Rm  and  are  provided  by  the  position  analysis  (Equations  A. 9  and  A, 23}, 
(d'Ti,,)/,  is  given  by  Equation  2,35  in  Step  3,  and  (Afi,,)*  is  provided  by  Equation  2.8  in 
Step  3  of  Section  2,1. 

Step  5r  From  the  position  analysis  (Equation  A. 26),  the  error  in  the  body-object 
displacement  vector  in  Earth  coordinates  is 

(Af),,),,  ^  {AfiA'  +  (2.41) 

Differentiating  (Ar^j)e  with  respect  to  time  yields  the  time  rate  of  change  of  the  error  in 
the  body-object  displacement  vector  (in  the  Earth  frame), 

=  (AF,,)„  (AFi,)„.  (2.42) 

The  full  expression  for  (fjrwjf  takes  the  form 

I  )r  =  R-l  If  +  +  6Rl  { A4 ),  +  Rli^Fu .  (2.43) 

Note,  the  first  two  terms  in  the  above  expression  mimic  the  result  for  (dFiJr  in  the  position 
analysis  (Equation  A. 27).  The  terms  /?,«,  (Afv,),,  and  (firi-i),  are  provided  by  the 

position  analysis  (Equations  A.1 1 ,  A.28,  A,8,  and  A.22  ),  is  given  by  Equation  2.39 
in  Step  4,  and  (An^,),  is  provided  by  Equation  2,11  in  Step  4  of  Section  2,1 ,  The  term 
Rf^  can  be  obtained  from  Equation  2,16  in  Step  5  of  Section  2,1 .  The  term  SR'j^  can  be 
derived  using  the  results  from  the  position  analysis  (Equation  A.28), 

6Ri,  =  SiltJi.tr  ■  (2.44) 
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Taking  the  time  derivative  of  this  expression  yields 


(2.45) 


The  term  is  provided  in  the  position  analysis  (Equation  A.29)  and  its  time  derivative 
is 


0  —6  0^2 

—  6&:^  0  S<l)\ 

662  D 


(2.46) 


Also  from  the  position  analysis  (Equation  A.30), 


■  0  ■ 

'  c. 

-  —So 

1 

-1-  <SA 

0 

.  . 

J 

0 

i 

.  . 

which  has  the  time  derivative, 


On.  Or, 

■  5'^  ' 
(1 

Stj) 

~  0  ■ 
1 

+  (iA 

■  ' 

0 

1 

1 _ 

I 

1 

J 

t 

0 

i 

-So 

(2.47) 


(2,48) 


Note,  the  last  two  terms  in  the  above  expression  mimic  the  result  for  {6^)1  in  the  position 
analysis  (Equation  A, 30). 

Step  6;  From  the  position  analysis  (Equation  A. 31 ),  the  error  in  the  object  position  vector 
in  Earth  coordinates  is 

(f)).  =  (r,).-(rf;v)..  (2.49) 

Differentiating  (7)^  with  respect  to  time  yields  the  time  rate  of  change  of  the  error  in  the 
object  velocity  vector  (in  the  Earth  frame), 


(AFJ,  =  (?,)- (2-50) 

which  is  equal  to 

^  (2-51) 

The  term  )r  is  known  fronn  Equation  2,43  in  Step  6  above.  The  term  (^Ti,)^  is  the  error 
in  the  vehicle  velocity  vector  expressed  in  the  Earth  frame,  and  must  be  determined. 
From  the  position  analysis  (Equation  A, 33),  we  know 


(4rt,)c  =  t  /j)4> 

■  ■ 

-■%s. 

T  (^’pEi'  T  /3)Cp(T\ 

■  -.5^  ■ 

Cx 

+  6h 

1 - 

1 _ 

^■'4  . 

U 

c 
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The  time  derivative  o1  yields 


dip 


l/'ni  +  h](Si^5xi\  —  C^Cx<p}  —  +  h)S^Cx 

-(r„^  +  +  C^)  -  [rn.  +  h)S^Sx 

“  ( ^  “1“  J^O  O  ~\~  (/  Hr  /i ) 


Hh 


(^ew  +  h)[—Ci^C),\  +  S^Sxip)  -  +  h}Cr^,Sx 

^  c  u;  +  ^0  {  ^■'  H"  Sf^C\ ^"^ )  -|-  [  f'c  Li.  "h  ^0  Co  (  A 

L) 


H- 


'  -  .v;,o  ■ 

'  -5^Ca  ■ 

6h 

•h  (rna  -1-  h)S(;) 

1 

e 

■  -Sa  ■ 
('■a 

"h  d// 

1 

O 

o  o 

E _ 

0 

e 

L  J 

(2.53) 


where 

r,,„.  =  5^CV^>(  1  -  (2.54) 

and 

r„,  =  :j(l  -  t^)r,M  -  (2.55) 

When  computing  the  velocity  error  sensitivities,  we  will  want  to  express  the  sensitivity  of 
the  object  velocity  error  to  north  and  east  position  (and  velocity)  errors.  Therefore,  it  is 
useful  to  express  in  terms  of  errors  in  north  and  east  coordinates.  Error  in  latitude 
(rt$>)  maps  to  error  in  north  position  as  to,v  =  S6{r„f  +  h).  Longitude  error  (^A)  maps 
to  east  position  error  as  =  C^iSXit'ac  +  When  these  substitutions  are  made,  the 
resulting  expression  for  (A’Ty,)^  is: 


.%CiA  -  C^Sf,o 

“h 

1 

1  1 
■Jj  1^1 

+  d/i 

■  -CXSj  -  5iCA0  ' 
tTiiC'AA  - 

c 

0 

F. 

■  -AVCa  ^ 

1 

1 

L _ 

C'a 

(} 

<• 

6V 

(2.56) 


2.3  Dynamic  Errors 
2.3.1  Tracker  errors 

The  object  velocity  error  incurred  as  a  result  of  the  tracking  algorithm  is  given  by 

(rlT),  4- (riVT),Ar,  (2.57) 

where  AT  is  the  true  time  latency  of  the  system.  The  term  (rtiT)f.  represents  any  error 
the  tracking  algorithm  adds  to  the  object  velocity  (e.g,.  smoothing  errors),  The  term  (rtfT)^ 
represents  the  error  in  the  acceleration  used  in  the  tracking  algorithm  to  extrapolate  the 
object  position.  These  errors  are  expressed  in  the  Earth  frame, 
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2.3.2  Tfming  errors 

The  object  velocity  error  incurred  as  a  resuit  of  a  timing  error  (^i)  is  given  by 


(2.58) 


where  (aj^  is  the  object  acceleration  expressed  in  the  Earth  frame.  Timing  errors  can 
include  time  stamp  error,  time  synchronization  errors  {e.g.,  host-to-network  error),  time 
latency  errors,  etc. 

2.4  Total  Aerospace  Object  Velocity  Error 

The  totai  error  in  the  object  velocity,  expressed  in  the  Earth  frame,  is  the  sum  of  the 
previous  results: 


- \r^ - 

dyciatnic 


The  total  error  can  be  partitioned  into  data  registration  errors  and  dynamic  errors  as 
indicated  above. 
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3  THE  VELOCITY  ERROR  SENSITIVITIES 

With  the  object  velocity  vector  completely  described  in  Section  2.1  and  the  errors 
derived  in  Sections  2. 2-2, 4.  the  velocity  error  sensitivities  can  be  derived.  The  ability  to 
express  the  sensitivities  explicitly  is  a  direct  result  of  the  linear  error  analysis  employed 
in  Section  2,  We  are  interested  in  defining  the  sensitivities  of  the  object  velocity  error  to 
the  variOLts  individual  error  sources  {i.e.,  position  errors,  tilt  errors,  time  synchronization 
errors,  etc.)  The  total  velocity  error  can  be  compactly  written  as 


(3,1) 


where  5  is  a  3  x  A’  sensitivity  matrix  (A'  is  the  number  of  individual  system  error  sources 
considered)  and  e  is  the  A'  x  1  system  error  vector.  This  form  is  advantageous  since  the 
covariance  of  the  track  velocity  error  can  be  directly  obtained  from  the  expression 


cov{6r,X  =  5cov{e)  S'^. 


(3-2) 


Similarly,  the  covariance  of  the  track  velocity  error  for  groups  of  error  contributors  (e.g., 
all  of  navigation)  can  be  easily  calculated  by  considering  the  appropriate  submatrices  of 
S.  The  elements  of  5  provide  information  on  "how  fast”  the  components  of  track  velocity 
error  change  as  the  individual  error  terms  change.  The  three  rows  of  S  correspond  to 
the  ,T-,  y-,  and  z-components  of  the  WGS  84  coordinate  system  and  each  column  of  S 
corresponds  to  an  element  of  the  error  vector  In  the  expressions  below,  the  notation  5^, 
represents  a  single  column  of  the  sensitivity  matrix  and  denotes  the  track  velocity  error 
sensitivity  to  the  error  e. 

In  what  follows,  the  north  and  east  velocities  of  the  vehicle  need  to  be  converted  to  time 
rates  of  change  of  latitude  and  longitude: 


(3,3) 


and 


(3,4) 


To  north  position  error: 


(3.5) 
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To  east  position  error: 


^  A 


« 


+ 


Hi 


^  c  uj  ^0 


^4^  ^  ^  hij  i 


Cp(j’cic  + 


^  I^J-) 


So^nt 


C^iXn^i 


+ 


+ 


J  t 


'  r  .■c-TiT^"'  '-■i  /  ■* 

f'_/"'.AJ_  A/vj  h-  r; .  ■^cu.’tJiA?' 


-C  +  A.^ 


{3.6) 


To  height  error: 


To  north  velocity  error: 


^Dh  = 


f ■'  *-^  X  A  ^  A  ^ 
C^CaA  -  S^S^<p 


- 

C^<i> 

-*  f 

'  -5, /A  ■ 
-5^Sa 

{] 

.  - 

r„s  +  h 

c 

(3.7) 


(3.8) 


To  east  velocity  error: 


-5a  ■ 
C’a 
0 


+ 


Hi 


6(!>Aj  (i;  ,^ 


(3.9) 


To  vertical  velocity  error: 


C^C\ 

C^Sx 
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To  navigation  tilt  (north)  error; 


5^,.,  =  Rijri; 


(J 

r  0  ■ 

+ 

b 

- 1 

<] 

(3.11) 


To  navigation  tilt  (east)  error; 


— 1 

[> 

c- 

— j 

D 

+ 

0 

b 

(3.12) 


To  navigation  tiit  (down)  error; 


bf-hj 

1 — 

1 _ 

b 

0 

(3.13) 


To  lever  arm  (6i)  error; 


■S’Ait  = 


- 


(3.14) 


To  lever  arm  ([i^}  error: 


_  i?T,jT 


(3.15) 


To  lever  arm  (&;,)  error: 


SSH,  = 


- 


(3.16) 
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To  sensor  misalignment  (vi ): 


Si. 


u 

1 

O 

_ 1 

+  Rl} 

ly 

To  sensor  misalignment  (fi'a): 


Sixjii  =  ^icSL 


0 

+  Rl.sl 

LI 

h 

JJ 

.  ^ 

_ I 

To  sensor  misalignment  (  ?■'«): 


^^'Sth2 

+  RI.RI, 

0 

(« 

0  J 

J  b 


(3.17) 


(3.18) 


(3.19) 


To  sensor  range  error; 

Sn/t  =  S.lfS.liRtts 


~Sfj, 


+  RifilH.,, 

Considering  only  range  and  range  rate  measurements 


^,Co,Co,  -  (hS,...S,, 


Sfif)  —  RffRhi  R-b-t 


t\C\ 

Sf), 


(3.20) 


(3.21) 


To  sensor  bearing  error: 


RlRlRb.{^ 

pG.CV.  -  pthS,yJ\  - 

L) 

A 

(1 

Considering  only  range  and  range  rate  measuremenls, 


■  ■ 

+  RlRiR,.p 

■  ■ 

Rl  n7,n,,.p 

t.) 

0 

(3.22) 


(3,23) 
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To  sensor  elevation  error: 


S —  Rjf.  R  p 

'  -Se,C^;' 

-Co, 

5 

— pC-ij,  -t- 

Considering  only  range  and  range  rate  measurements, 


(3.24) 


1 

L _ 

ff 

a.. 

{3.25) 


To  sensor  range  rate  error: 


C^->Co, 

So. 


S^-.C 


To  sensor  bearing  rate  error: 


Rb,p 


Sip,Ci)^ 

0 


To  sensor  elevation  rate  error; 


-Q. 


To  tracker  position  errors: 


0 

0 

n 


(3.2G) 


(3.27) 


(3.28) 


{3.29) 
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To  tracker  z-veiocity  error: 


(3-30) 


To  tracker  i^-vejocity  error: 


(3.31) 


To  tracker  ^-velocity  error: 


(3.32) 


To  tracker  .r-acceleration  error: 


5, 


daj; 


{3.33} 


To  tracker  u-acceleration  error; 


3  Ha,.  — 


AT 


(3.34} 


To  tracker  ^^acceleration  error: 


rlrj.- 


AT 


(3.35) 


To  timing  errors; 


a- 


(3.36) 
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4  RESULTS 

4.1  Scenario  Description 


To  develop  some  physical  intuition,  we  will  consider  a  simple  tracking  case  to  identify 
the  important  sensitivities  and  quantify  the  relative  importance  of  the  error  contributors. 
These  results  can  be  compared  to  the  results  for  position  error  for  additional  insight.  The 
scenario  {see  Figure  4.1)  consists  of  a  sea-based  sensor  tracking  an  inbound  cruise 
missile,  The  operational  parameters  for  the  scenario  are  listed  in  Table  4,1 ,  The  sensor 
is  mounted  to  a  ship  that  is  traveling  at  1 0  kts  on  a  NW  course  off  the  coast  of  South 
Korea.  The  sensor  frame  is  directionally  aligned  with  the  vehicle  (body)  frame.  The  cruise 
missile  target  is  traveling  at  600  kts  due  East  and,  at  the  particular  instant  considered,  is 
undergoing  a  1  -g  acceleration  in  the  East  direction. 

4.2  The  Sensitivities 

The  transpose  of  the  sensitivity  matrix  5  is  shown  in  Table  4.2.  These  sensitivities  were 
computed  using  the  simplifying  assumptions  that  the  sensor  frame  is  fixed  relative  to 
the  vehicle  body  frame  and  that  the  vehicle  tilt  (attitude)  rates  are  sufficiently  small.  In 
addition,  neither  bearing  rate  nor  elevation  rate  measurements  are  considered.  The 
numerical  entries  in  the  table  represent  the  amount  of  track  velocity  error  in  each  of  the 
three  axes  (xv,  y,,  of  the  WGS  84  coordinate  frame  that  are  generated  by  a  unit  error 
term.  For  example,  if  is  the  first  dimension  of  track  velocity  and  vehicle  east  velocity 
is  the  fifth  error,  then  a  value  of  -0,777  for  the  (1 ,5)  element  of  the  sensitivity  matrix 


NOT  TO  SCALE 


Figure  4. 1 :  Example  Tracking  Scenario. 
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Item 

Parameter 

Value 

Units 

Latitude 

deg  N 

36 

Longitude 

deg  E 

130.05 

Altitude 

m 

0 

North  velocity 

kts 

7.1 

East  velocity 

kts 

-7/1 

Vehfcle 

Down  velocity 

m/s 

0 

cpb  (roll) 

deg 

0 

(pitch) 

deg 

0 

V-b  (yaw) 

deg 

315 

Lever  arm,  (A/|,Jt, 

m 

-30 

Lever  arm. 

m 

Lever  arm,  ( 

m 

-2 

tp,,  (roll) 

deg 

0 

e,  (pitch) 

deg 

0 

■Ip,  (yaw) 

deg 

0 

Sensor 

Range  to  target 

nm 

50 

Range  rate  i 

m/s 

-308 

Bearing  to  target 

deg 

317 

Elevation  to  target 

deg 

2 

Network 

Latency 

sec 

5 

■Ts- velocity 

m/s 

-240 

y, -velocity 

m/s 

-194 

Aerospace 

-velocity  [ 

m/s 

0 

Object 

acceleration 

Q 

-0.78 

yg-acce]  oration 

9 

-0,63 

jjp-acceferation 

9 

0  1 

Table  4,1:  Scenario  Operational  Parameters 


means  that  every  1  m/s  of  east  velocity  error  causes  0.777  m/s  of  track  velocity  error  in 
the  negative  Xe -direction.  If  the  east  veJocity  error  is  1  kt,  then  its  contribution  to  track 
velocity  error  is  0.4  m/s  in  the  negative  ait, -direction. 

Insight  into  how  individual  errors  contribute  to  track  velocity  error  can  be  gained  by 
examining  the  root-sum-square  (RSS)  of  the  elements  of  the  sensitivity  matrix.  These 
values,  referred  to  as  “Unit  Radial  Track  Velocity  Error,"  are  provided  in  the  last  column 
of  Table  4.2.  The  values  of  unity  for  Unit  Radial  Track  Velocity  Error  corresponding  to  the 
three  vehicle  velocity  error  components  imply  that  1  m/s  of  vehicle  velocity  error  in  the 
north  or  east  or  down  direction  translates  one-for-one  into  1  m/s  of  track  velocity  error. 
This  mapping  is  independent  of  operational  conditions  and  scenario,  e.g.,  aerospace 
object  conditions.  The  same  holds  for  the  contributions  of  sensor  range  rate  error  and 
tracker  velocity  errors. 
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Error  Term 

Unit 

Sensitivjties  (m/s/unit) 

Unit  Radial 

WGS  84  Earth  Fixed  Axes 

Tracit 

Velocity 

■^c 

Error 

North  position 

m 

0.000 

0.000 

0.000 

0.0[)(1 

East  position 

m 

u.ono 

0.000 

0.000 

o.oou 

Height 

m 

U.OULf 

0.000 

Q.OUO 

0.000 

North  velocity 

m/s 

0J7<]  ■ 

-0.450 

1  n.so!) 

1.000 

East  veJocity 

m/s 

-0.777 

-0.650 

O.OOff 

1.000 

Down  velocity 

m/s 

-0,521 

0.619 

0.588 

1.000 

North  tilt 

mrad 

U.113 

-n.122 

-0.150 

0.211 

East  tilt 

mrad 

0.126 

-0.158 

-0.126 

0.225 

Down  tilt 

mfad 

-0.109 

0T45  , 

-‘0.240 

0.508 

Lever  arm,  bi 

m 

11.000 

0.000 

0.000 

Lever  arm, 

m 

o.oon 

0.000 

0.000 

O.OOU  ! 

Lever  arm,  63 

m 

0.000 

0.000 

U.OOO 

o.ono 

Sensor  IVlisafigrtmenL 

mrad 

0.122 

DT5U 

IT211 

Sensor  Misalignment,  th 

mrad 

-0.126 

0.158 

0.126 

il.225 

Sensor  Misalfgnment; 

mrad 

0.109 

-0.145 

0.249 

0.308 

Range  to  object 

m 

0.000 

0.000 

0.000 

r  (j.noo 

Range  rate 

m/s 

0-759 

0.050 

0.047 

i.wo 

Bearing  to  object 

mrad 

■0.109 

0.145 

-0.249 

(J.8O8 

Elevation  to  object 

mrad 

0,109 

-0.134 

-0.181 

U.30S 

Tracker  avposition 

m 

U.UOO 

0.000 

0.000 

O-OOtl 

Tracker  7/4, -position 

m 

U.OOO 

0.0(J() 

o.noo 

o.oon 

Tracker  .if- position 

m  1 

0.000 

0.000 

0.000 

O-Oan 

Tracker  .3.:^“  veJodty 

m/s 

i.oon 

o.oou 

0.000 

1.000 

Tracker  T/ft- velocity 

m/s 

o.uou 

i.ann 

O.OUO' 

L.rjoo 

Tracker  ^t--'^©J<5city 

m/s 

o.noo 

O-OUO 

LOOO 

i.non 

Tracker  -acceleration 

m/s^ 

5.DO0 

o.oorj 

o.uou 

5.000 

Tracker  73^ -acceleration 

m/s^ 

ITODO 

5.000 

0.000 

5.000 

Tracker  -acceleration  | 

m/s^ 

o.ono 

o.uou 

5.000 

5.000 

Timing®  | 

sec 

^7.029  ; 

-ti.l84 

n.ooo 

Sl.821 

^Represents  the  sensflivity  of  any  time  error  such  as  time  synchronization,  time  stamp,  etc. 

Table  4.2:  Velocity  Error  Sensitivity  Matrix  (Transposed)  and  Unit  Radial  Track  Velocity 
Error 


Track  veJocity  error  is  not  sensitive  to  position-type  errors  such  as  navigation  position 
error,  lever  arm  errors,  sensor-measured  range  error,  and  tracker  position  errors. 

As  already  noted,  track  velocity  error  has  a  1:1  sensitivity  to  velocity-type  errors. 
Orientation-,  or  angular-type,  errors  contribute  a  small  amount  to  track  velocity  error 
for  this  scenario.  These  include  navigation  tilt  errors,  sensor  misatignments,  and 
sensor-measured  bearing  and  elevation  errors.  The  sensitivity  to  navigation  tilt  errors  is 
practically  identical  to  the  sensitivity  to  sensor  misalignments  due  to  the  small  lever  arm 
vector.  Likewise,  the  bearing  and  elevation  sensitivities  are  effectively  identical  because 
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Error  Term 

Unit 

Sensitivities  (m/$/unit) 

Unit  Riadj^l 

Track 

Error 

WGS  84  Earth  Fixed  Axes 

North  position 

m 

-0.450 

0.809 

I.OOO 

East  position 

m 

-0.777 

-0.680 

UJJOO 

1 .000 

Height 

m 

-0-521 

0.610 

0.58S 

L.ono 

North  velocity 

m/s 

0.000 

u.ooo 

G.OOO 

0.000 

East  velocUy 

m/s 

0.000 

0.000 

O-OOO 

o.ono 

Down  velocity 

m/s 

0.0  no 

0.000 

fl.OOO 

n.ooo 

North  tilt 

mrad 

-33.815 

86.670 

39.038 

68.388 

East  tilt 

mrad 

-87.720 

41,826 

87J92 

67.524 

Down  tilt 

mrad 

82J>8tl 

-4.8.40fl 

74.769 

92,468 

Lever  arm,  by 

m 

0.800 

0.137 

0.572 

1.000 

Lever  arm,  id 

m 

■0.27-1 

-0T78 

0.572 

1,000 

Lever  arm, 

m 

0.521 

-Dim 

-0.588 

1.000 

Sensor  MisalignmenL 

mrad 

33.S12 

-86.67-5 

’-80.027 

63.33G 

Sensor  Misalignment,  i^'2 

mrad 

37.788 

-41.845 

-37.814 

67.558 

Sensor  Misalignment 

mrad 

-32.676 

43.514 

-74.785 

92.488 

Range  to  object 

m 

U.fi-'iO 

0.047 

1.000 

Range  rate 

m/s 

d-OOfj 

O.MO 

0,000 

0.000 

Bearing  to  object 

mrad 

-43.5M 

74.785 

92.488 

Elevation  to  objecl 

mrad 

-  Tid.fiDr 

55.247 

54.279 

92.544 

Tracker  iTTf^-posilEcn 

m 

t.ouo 

0.000 

0.000 

1.000 

Tracker  position 

m 

[kOUO 

i.OOO 

0,000 

1 .000 

Tracker  s-^-position 

m 

o.ouo 

0,000 

LODO 

1.000 

Tracker  i^.-velocity 

m/s 

5.000 

0,000 

0.000 

r».o(){i 

Tracker  -'velocity 

m/s 

Q.OOO 

5.000 

0.000 

5.0G0 

Tracker  ‘Velocity 

m/s 

0.000 

0.000 

5.000 

5-000 

Tracker  2:e-acceleration 

m/s^ 

I2.50U 

0.000 

o.ono 

12.500 

Tracker  t/c -acceleration 

m/s^ 

u.uoo 

12.500 

04)00 

12,500 

Tracker  5:^ -acceleration 

m/s^ 

o.ono 

0.000 

12.500 

12.500 

Timing^ 

sec 

-‘^77,^37 

-225.279 

-0.005 

357-770 

^Represents  the  sensitivity  of  ary  time  error  such  as  Itme  synchronisation,  time  stamp,  etc. 


Table  4.3:  Position  Error  Sensitivity  Matrix  (Transposed)  and  Unit  Radial  Track  Error 


of  the  low  elevation  of  the  target.  The  sensitivity  to  timing  errors  dominates  the  sensitivity 
matrix.  This  component  is  only  dependent  upon  the  acceleration  of  the  aerospace  object. 
The  greater  the  acceleration,  the  stronger  the  sensitivity. 

For  comparison,  the  track  position  error  sensitivity  matrix  is  shown  in  Table  4.3.  The 
first  thing  to  note  is  that  track  position  error  is  generally  much  more  sensitive  to  the 
various  error  terms.  For  this  scenario,  angular-type  errors  and  timing  errors  dominate 
the  position  error  sensitivity  matrix.  There  are  some  similarities  between  the  velocity 
error  sensitivities  and  the  position  error  sensitivities.  Specifically,  we  recognize  that 
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the  position  error  sensitivities  to  navigation  position  errors  and  sensor- measured  range 
errors  are  equivaJent  to  the  velocity  error  sensitivities  to  navigation  veiocity  errors  and 
sensor- measured  range  rate  errors,  which  are  just  the  time  rates  of  change  of  navigation 
position  and  range. 


4,3  System  Error  Values 

A  nominal  error  budget  for  this  scenario  is  presented  in  Tabie  4.4^  The  errors  are 
assumed  to  be  normally  distributed  with  zero  mean.  For  this  analysis,  alf  errors  are 
assumed  to  be  independent.  However,  the  mefhodoiogy  presented  in  Section  2  does 
not  require  this  lo  be  true.  The  uncorreiated  model  was  chosen  for  convenience.  A 
detaiied  analysis  of  the  individuai  subsystems  is  necessary  to  develop  the  more  generai 
correlated  error  model. 


4,4  Track  Velocity  Error  Budget 

Combining  the  sensitivity  matrix  in  Tabie  4.2  with  the  nominai  error  budget  jn  Tabie  4.4 
leads  to  the  nominai  track  veiocity  accuracy  error  budget  dispiayed  in  Figure  4.2.  The 
bars  in  the  figure  are  color-coded  to  heip  visuaiize  the  roli^up  of  individuai  errors 
Navigation  errors  are  colored  in  light  blue,  sensor  errors  are  medium  blue,  tracker  errors 
are  green,  and  timing  errors  are  ,  Above  the  horizontai  dotted  line  are  the  individual 
error  contributors  and  beiow  the  line  are  the  groupings  of  error  contributors,  which  are 
coior-coded  accordingiy.  In  addition,  the  data  registration  error  rofi-up  (orange),  the 
dynamic  error  rolhup  (red),  and  the  totai  system  error  (black)  are  shown. 

The  track  veiocity  accuracy  is  quantified  in  terms  of  the  95%  radius,  which  represents  the 
radius  of  a  sphere  that  encompasses  95%  of  the  track  veiocity  error  distribution.  Unlike 
the  case  of  track  position  error,  this  95%  sphere  does  not  have  a  physical  interpretation 
since  the  radius’  is  in  units  of  m/s.  Referring  to  Figure  4.2,  the  iength  of  the  bar  for  each 
error  term  iisted  in  the  chart  reflects  the  contribution  to  track  velocity  error,  expressed  as 
95%  radius,  if  that  error  term  was  the  only  error  In  the  system.  For  instance,  if  range  rate 
were  the  only  error,  the  95%  radius  for  the  system  would  be  about  20  m/s. 

This  chart  aiso  lilustrates  the  cumulative  effects  of  several  errors.  For  exampie^  the 
sensor  range,  range  rate,  bearing,  and  elevation  error  contributors  combine,  or  roll  up, 
to  give  a  total  MEASUREMENT  error  contribution  of  approximately  20  m/s  95%  radius, 
which  is  dominated  by  the  range  rate  error.  Simiiarfy,  the  tracker  position,  velocity,  and 
accelemtion  error  contributors  combine  to  yieid  a  total  TRACKER  error  contribution  of 
approximately  50  m/s  95%  radius.  The  total  system  error  of  53.6  m/$  95%  radius  is 
composed  of  a  DATA  REGSTRN  (NAVIGATION  -f  LEVER  ARM  +  SENSOR)  contribution 
of  19.4  m/s  95%  radius  and  a  DYNAMIC  (TRACKER  +  TIME)  contribution  of  50.9  m/s 
95%  radius.  For  comparison,  the  totai  system  error  for  track  position  is  781.3  m  95% 
radius  for  this  scenario.  Consequentiy,  it  wouid  take  aimost  15  seconds  for  this  velocity 
error  (53,6  m/s)  to  grow  into  an  equivalent  position  error  of  781  m.  This  relatively 
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Item 

Error  Parameter 

Unit 

1 -Sigma  Value 

North  position 

m 

200 

East  position  , 

m 

200 

Height  i 

m 

20 

North  velocity 

m/s 

0.25 

East  velocity 

m/s 

0.25 

Vehicle 

Down  velocity 

m/s 

0.10 

North  tiit 

mrad 

0,50 

East  tiit 

mrad 

0,50 

Down  tilt 

mrad 

0,50 

Tever  arm,  (ftrL)i>„ 

m 

0.10 

Lever  arm, 

m 

0.10 

Lever  arm, 

m 

0-10 

Sensor  Misalignment,  &'i 

mracH 

1 

Sensor  Misalignment,  % 

mrad 

1 

Sensor  Misalignment,  ^'>3 

mrad 

1 

Sensor 

Range  to  object 

m 

100 

Range  rate 

m/s 

10 

Bearing  to  object 

mrad 

2 

Elevation  to  object 

mrad 

2 

Tracker  .r,, -position 

m 

10 

Tracker  t/r-position 

m 

10 

Tracker  ir^-posltion 

m 

10 

^Tracker  ,-<,y -velocity 

m/s 

15 

Tracker 

Tracker  -velocity 

m/s 

15 

Tracker  Sf-veiocity 

m/s 

15 

Tracker  ;c,>-acceleration 

m/s^ 

2 

Tracker  ^^-acceleration 

m/s^ 

2 

Tracker  Sc-acceleration 

m/s^ 

2 

Host-to-host 

sec 

0.25  ' 

Host-to-network 

sec 

0.25 

Timing 

Time  stamp 

sec 

0.25 

Time  latency 

sec 

0.25 

Table  4.4:  Nominal  Scenario  Error  Budget 
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N  position 
ii  position; 

Heighr 
H  velcicit>'  I 
E  vektity  I 
D  veloGTt)-  p 
North  fc1  Y 
Efti«  tin  [ 
Down  SilL  p 
Laver b1  p 
(.Q^-ar  b2 
Laver b3 
Sensor  psi1 
SePEor  sMffi 
Sansor  psa 
.Range . 
Harm  a  r.3^ 
££anng 
Elavatibn 
Trar^sr  ji.-ppE 
Trader  y-  pos 
T  racJier  i-paa 
T  racka:  K-vel 
Tretkai  y-vel 
Tracker  7. 1^1 
Traclflsf  x-acc 
Tracker  y-acE 
Tracker  z-ace 
Hosl-lo-host 
Hoet-to-net 
Tima  stamp 
Time  talancy  _ 
PQS5TJON  - 
VELOCITY 
T^LT 

NAVIGATION 
L^VER  ARM 
.MEASUREMENT 
MISALIGNMENT 
SENSOR 
1  HACKER 
TIME 

OATA  REGSTRW 
DVMAMIC 
TOTAL 


1^5%  radius  (rn/s) 


Figure  4.2:  Nominal  Track  Velocity  Accuracy  Error  Budget. 


long  time  implies  that  track  velocity  error  is  not  as  important  as  track  position  error  for 
successful  aerospace  object  tracking. 

The  error  budget  results  provide  an  easy  way  to  identify  the  major  contributors  to  total 
track  velocity  error,  at  both  the  subsystem  level  and  at  the  individual  error  component 
level.  For  this  scenario,  dynamic  errors  dominate  the  total  system  error.  These  dynamic 
errors  are  composed  of  tracker  and  timing  error  contributions.  Specifically,  the  tracker 
velocity  and  acceleration  error  contributions  are  the  primary  cause  for  the  large  tracker 
error  contribution.  Therefore,  improvements  to  the  tracker  algorithms  would  yield  the 
greatest  benefit  in  terms  of  reducing  track  velocity  error.  Additional  benefit  could  be 
gained  with  an  improved  sensor  measurement  of  range  rate,  improvements  to  the  vehicle 
navigation  system  would  have  little  effect  on  track  velocity  error  as  evidenced  by  their 
negligible  contributions  in  the  chart. 
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5  CONCLUSIONS 


This  Technical  Report  has  provided  the  details  of  an  analytical  technique  to  understand 
and  quantify  how  errors  in  navigation,  sensor,  timing,  and  tracker  affect  the  accuracy  of 
track  velocity  estimation.  The  analysis  approach  derives  the  equation  for  the  velocity 
of  an  aerospace  object  in  the  WGS  84  coordinate  frame  and  determines  how  data 
registration  and  dynamic  errors  affect  the  estimate  of  track  velocity.  An  examination 
of  the  components  making  up  the  track  velocity  vector  reveals  three  groups  of  error 
contributors: 

•  Vehicle  position,  velocity,  and  attitude 

•  Sensor  measurements  and  misalignments 

•  Timing  (e.r?.,  time  stamp,  time  synchronization,  time  latency,  etc.) 

•  Tracker 

Comparing  these  results  for  velocity  error  to  those  for  position  error,  we  identify  new 
contributors  to  track  velocity  error: 

•  Vehicle  velocity 

•  Sensor  rate  measurements  (if  available} 

Vehicle  tilt  rate  {attitude  rate)  errors  do  not  enter  into  the  track  velocity  error  in  this 
analysis  because  of  the  simplifying  assumptions  made  in  Section  2.  After  the  error 
mechanisms  were  identified,  a  first-order  linear  perturbation  analysis  was  employed  to 
develop  the  relationship  between  track  velocity  error  and  the  errors  identified  above. 

A  key  result  of  the  analysis  is  the  determination  of  the  sensitivities  of  track  velocity 
accuracy  to  the  navigation,  sensor,  timing,  and  tracker  errors.  With  the  sensitivities 
established,  a  nominal  tracking  scenario  was  considered  and  a  detailed  track  velocity 
error  budget  generated.  The  error  budget  results  provided  an  easy  way  to  identify  the 
major  contributors  to  total  track  velocity  error,  at  both  the  subsystem  and  individual 
error  component  level.  This  insight  shows  where  improvements  in  the  system  could  be 
made  in  order  to  yield  the  greatest  benefit  in  terms  of  track  velocity  error.  The  results 
further  indicate  that  track  velocity  error  contributors  are  less  important  than  track  position 
error  contributors  for  successful  aerospace  object  tracking.  As  a  result,  the  Integrated 
Architecture  Behavior  Model  (lABM)  development  activities  should  focus  on  capturing 
the  position  error  contributors  for  the  Configuration  05  build, 

Finally,  the  results  from  this  analysis  can  be  used  to  investigate  track  vefocity  error  and 
will  be  incorporated  into  the  JSSEO  Parametric  Track  Accuracy  Model  thereby  enhancing 
its  scope  and  capability.  For  example,  we  are  using  the  error  model  and  the  results  of 
the  analysis  runs  to  develop  test  cases  to  assist  with  the  tABM  verification  and  validation 
effort.  The  error  models  provide  the  mechanism  to  introduce  the  correct  level  of  bias 
into  each  system,  and  the  results  of  the  error  model  analysis  provide  examples  of  track 
velocity  error  that  the  lABM  can  be  validated  against. 
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APPENDIX  A:  SUMMARY  OF  POSITION  ERROR  ANALYSIS  RESULTS 
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A.1  The  Earth 


We  consider  the  WGS  84  ellipsoidal  model  of  the  Earth  geoid  (see  Table  B.1}.  The 
geodetic  coordinates  of  a  point  B  with  respect  to  this  datum  will  be  denoted  as  {o.  A,  /?), 
which  are  geodetic  latitude,  longitude,  and  altitude  of  point  B  above  the  ellipsoid.  The 
Conventional  Terrestrial  System  (CTS)  associated  with  the  WGS  84  datum  will  form  the 
right-handed  Cartesian  Earth  frame,  or  <;-frame. 


A.2  The  Vehicle 


We  consider  a  vehicle,  fixed  or  moving,  in  the  vicinity  of  the  Earth,  We  define  a  point  B, 
fixed  to  the  vehicie,  representing  the  center  of  the  vehicle’s  Inertial  Navigation  System 
(INS)  with  geodetic  coordinates  {(p,  A,  /f)  as  above.  The  vector  from  the  center  of  the 
Earth  to  point  B  on  the  vehicle,  expressed  in  the  r-frame,  is  given  by 


(r,,„ 

(■r,„  -h  h}C^S,,, 

[(1  -  -h  /i]  J  ^ 


(A.t) 


where  C’,.  and  are  shorthand  for  cosri  and  sino,  and  rp„, 
curvature  given  by 


is  the  east-west  radius  of 

(A.2) 


A.3  The  Aerospace  Object  Position  Vector  Processing  Chain 

The  following  is  a  summary  of  the  nominal  (error-free)  processing  chain  for  the  aerospace 
object  position  vector.  The  starting  point  are  the  sensor  measurements  of  the  object  and 
the  ending  point  is  the  object  position  vector,  which  is  the  vector  from  the  center  of  the 
Earth  to  the  aerospace  object,  expressed  in  Earth  coordinates. 


Step  1 :  Assume  that  the  sensor  system  associated  with  the  vehicle  measures  the  range 
p,  the  bearing  and  the  elevation  0,  of  the  object  with  respect  to  the  sensor  frame 
(,s-frame).  This  provides  the  sensor-object  displacement  vector  in  sensor  coordinates: 


{Ar\,i)s  =  p 


{A.3) 


Step  2:  The  sensor-object  displacement  vector  is  transformed  to  body  coordinates: 

(A. 4) 

The  sensor-body  rotation  matrix,  /4,,  is  given  by 

%  =  (A.  5) 
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where  are  the  three  EuJer  angles  describing  the  rotation  from  the  sensor 

frame  to  the  body  frame,  and  and  are  the  standard  trio  of  single-axis  rotation 

matrices  defined  as 


= 


1 

il 

0 

0 

slnB 

G 

—  sin 

.8  con 

cosO  0 

-  sill  (I 

Q  1 

0 

sin  (?  0 

ros  0 

t: 

i9 

SlLL(f  0 

— 

sin^ 

0 

D 

IJ  1 

(A.6) 


Step  3:  The  body-sensor  lever  arm,  is  added  to  produce  the  body-object 

displacement  vector  in  body  coordinates: 


(A.7) 


Step  4.  The  body— object  displacement  vector  is  transformed  to  local*level  coordinates; 

=  /?J(A4,)e,.  (A.8) 

The  local-ievel-body  rotation  matrix.  B,^,,  is  given  by 

Rm  =  (A.9) 

where  (i),,,  8^,  dj)  are  the  three  Euler  angles  describing  the  rotation  from  the  local-level 
frame  to  the  body  frame. 

Step  5:  The  body-object  displacement  vector  is  transformed  to  Earth  coordinates: 

(Afh,),.  = (A,  10) 
The  Earth-local-level  rotation  matrix,  R.u,  is  given  by 


B, 


/  7T  \ 

-,Va  C,,  “ 

i-r.  -0}  IMX\=. 

r;,  G 

-{'oCa  -Tr^5A  -.9^ 

(A, 11) 


Step  6;  The  navigation  system  position  vector,  (rj,),.,  is  added  to  produce  the  object 
position  vector  in  WGS  84  coordinates: 


(rj),  -  (n,),.  -I-  (4)r- 


(A.12) 
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A, 4  The  Position  Vector  Chain  o1  Errors 


The  same  steps  as  in  the  preceding  section  are  used  to  derive  Ihe  measured  quantities 
from  the  fundamental  error-free  quantities.  Measured  quantities  will  be  denoted  with 
and  error  quantities  wiii  be  denoted  with  a  preceding  6. 

Step  1:  Error  in  the  sensor-object  dispiacement  vector  in  sensor  coordinates: 

(A.  13) 


■  .5^aX'v.  ■ 

-\-dp 

- 

D 

_ 

> 

-S,.  J 

Step  2:  Error  in  the  sensor-object  dispiacement  vector  in  body  coordinates: 


—  6  ■ 


0  Sc- 

3 

4^1 

~  (.V I  ( 1 

-  ‘V's  ^ 

(A.15) 

(A.  16) 
{A.17) 

{A.18) 


Step  3:  Error  in  the  body-object  dispiacement  vector  in  body  coordinates: 

^  {A.  19) 

hj  —  A'f'bH  jH ( A . 20) 


Step  4:  Error  in  the  body-object  displacement  vector  in  iocaMevei  coordinates: 


mb  = 


{)Q{,t 


f)  -^60^. 

IJ  6$s 

— dW.v  0 


■  I  ■ 

0 

-Sh^ 

<yh: 

™  ±\ 

0 

-i-  df^y, 

+  <SL'|j 

CA 

h 

0 

b 

() 

(A.21) 

(A.22) 

(A.23) 

{A.24) 

(A.25) 


A-4  7.2.6.4-TR(4-001  ).1 ,0  JSSEO-04041 3 


step  5:  Error  in  the  body-object  displacement  vector  in  Earth  coordinates: 


^  Jt!  ‘ 


Si'h,  = 


0  rips  -Sft2 
-(^03  0  (fe| 

(bj  -fei  0 


{dS)t  - 

= 

^  0  ■ 
1 

“h  6X 

■  ■ 

(1 

. 

j 

0 

l 

.  . 

Step  6:  Error  in  the  aerospace  object  position  vector  in  Earth  coordinates: 

(ri)e  =  (''"Vlf  + 

r 

(f> '>.),:  =  (f,.,  + 

where  r,„  is  the  north-south  radius  of  curvature  given  by 

r„,  =  u,(l  -  e-)(l  -  -  e^)fi  - 


■  -s^c\ ' 

'  -Sx  ■ 

■  C^Cx  ' 

~S^s>Sx 

+  (l^ra-  + 

a 

+  (T?P. 

C^Sx 

CU 

e 

0 

4? 

(A.26) 

{A.27) 

{A.28) 

(A.29) 

(A.30) 


(A, 31) 
(A.32) 

{A.33) 

(A.34) 
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APPENDIX  B:  WGS  84  ELLIPSOID  PARAMETERS 
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Consider  an  ellipsoid  model  of  Ihe  Earth  geoid,  in  particular,  the  one  used  as  the 
geodetic  datum  by  the  World  Geodetic  System  1 984  (WGS  84).  Parameters  related  to 
this  ellipsoid  are  listed  in  Table  B.1. 


Parameter 

Notation 

Formula 

Va(uQ 

Angular  velocity 

7.2921 1 51467  x  10'^  rad/sec 
15.04106718  deg/hr 

Semimajor  axis 

6378137  m 

Flattening 

/ 

1/298.2572236 

Semiminor  axis 

h. 

«.(!  -/} 

6356752.3  m 

Eccentricity 

e 

^/7(2^ 

0.08181919084 

Eccentricity  squared 

0-00669437990 

Axis  ratio 

0.9966471 893 

Axis  ratio  squared 

l-e- 

0.9933056199 

Linear  eccentricity 

E 

ea^ 

521654.01  m 

Minor  eccentricity 

e' 

e/(l-/) 

0.08209443796 

Table  B.1:  WGS  84  ellipsoid  parameters 


B  2 
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